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Abstract — The Mars Science Laboratory entry capsule is used 
as an example to demonstrate how a blunt body of revolution 
must be treated as asymmetric in some respects when flying 
at a non-zero trim angle of attack. A brief description of the 
axisymmetric moment equations are provided before solving a 
system of equations describing the lateral-directional moment 
equations for a blunt body trimming at an angle of attack. Sim- 
plifying assumptions are made which allow the solution to the 
equations to be rearranged to relate the roll and yaw stability 
with sideslip angle to the frequency of oscillation of the vehicle 
body rates. The equations show that for a blunt body the roll 
and yaw rates are in phase and proportional to each other. The 
ratio of the rates is determined by the static stability coefficients 
and mass properties about those axes. A trajectory simulation is 
used to validate the static yaw stability parameter identification 
equation and a simple method of identifying the oscillation 
frequency from the body rates. The approach is shown to 
successfully extract the modeled yaw stability coefficient along 
a simulated Mars entry. Mars Science Laboratory flight data 
results are presented from earlier work which indicate that 
results from both the validation case and flight data are in 
agreement with preflight predictions. A brief discussion of 
the dynamic stability is also provided. Trimming at a non- 
zero angle suggests that the typical axisymmetric models of 
the dynamic stability coefficients should be modified. However, 
further experimental or computational work must be done to 
separate damping due to body rates and wind relative rates 
before the correct lifting formulation would affect simulation 
results. 
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1. Nomenclature 


Symbols 


A 

Sinusoid amplitude 

a, b, c 

Functional form coefficients 

C A 

Axial force coefficient 

Ca 0 

Derivative of Ca with respect to f3 

c La 

Lift curve slope with respect to a 


G h 

Roll stability with respect to j3 

Cm a 

Pitch moment slope with respect to a 

C N 

Side force coefficient 

Gn 0 

Derivative of side force with respect to f3 

Cn, 

Yaw moment slope with respect to (3 

c nr - c n . 

Yaw damping coefficient 

Cn r 

Yaw damping about body axis 

Cy 

Side force coefficient 

d 

Reference diameter 

Ixx ■ lyy 7 b zz 

Moments of inertia 

Ixy 7 lyz 7 z 

Cross products of inertia 

m , niEV 

Vehicle mass 

N damping 

Yaw damping moment 

p,q7f 

Body rates 

Po7fo 

Initial body rates (roll and yaw) 

Posc^osc 

Oscillatory components of p and r 

Qoo 

Freestream dynamic pressure 

r c 

Constant roll component 

fitted 

Sinusoid model fit to roll rate 
Reference area 

t 

Time 


Freestream velocity 

V X ,Vy,V Z 

Velocity components 

Vaci Zac 

Aerodynamic coordinate axes 

%cg i Ucg ) Z C g 

Center of gravity location 

% nose i Unose ? Z nose 

Nose location 

Greek 

a 

Angle of attack 

OL 0 

Initial angle of attack 

oi. r p 

Total angle of attack 

OtTrim 

Trim angle of attack 

P 

Angle of sideslip 

flo 

Initial angle of sideslip 

f3 Trim 

Trim angle of sideslip 

8 

Phase angle 

Poo 

Freestream density 

0 

Bank angle 

UJ 

Oscillation frequency 

W q 

Pitch oscillation frequency 

L O r 

Yaw oscillation frequency 


2. Introduction 

For planetary and Earth entry, blunt bodies of revolution have 
long been the geometry of choice for aeroshells used to slow 
payloads from hypersonic entry velocities to low supersonic 
or subsonic speeds where deceleration is augmented with 
parachutes and propulsion. Aeroshell heatshield designs 
have typically been made up of some combination of large- 
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angle cone and blunt sphere sections. The forebody shape 
is typically driven by a combination of heating and stability 
requirements depending on entry velocity and the planet 
atmosphere profile. Backshell geometries are typically ax- 
isymmetric as well with the specific shape frequently driven 
by payload constraints. Often, these bodies of revolution 
are flown with the center of gravity (eg) aligned with the 
axis of symmetry and therefore trim at zero degrees total 
angle of attack. This symmetry simplifies how aerodynamic 
forces and moments are modeled in any flight simulations 
used for mission design and landing site targeting. The 
static forces and moments can be tabulated as functions of 
total angle of attack only and decomposed into six-degree- 
of-freedom coefficients (with angles of attack and sideslip 
as the independent variables) for use in simulation. Often 
the dynamic stability derivatives are modeled as functions 
of total angle of attack as well. In contrast, a conventional 
airplane requires that aerodynamic coefficients be described 
independently about and along the roll, pitch and yaw axes 
(or other convenient six degree-of-freedom coordinate sys- 
tem). When an axisymmetric vehicle uses a radial center 
of gravity offset to fly at an angle of attack, relative to the 
approaching wind, axisymmetry is broken and in some ways 
the aerodynamic coefficients must be described more like an 
airplane. 

On August 5, 2012, the Mars Science Laboratory (MSL) 
entry vehicle successfully slowed the Curiosity rover during 
its entry, descent and landing in Gale crater on Mars. Both 
the aeroshell and entry architecture for MSL were derived 
from previous successful configurations used to land on Mars 
which have all been based on the Viking missions of the 
1970’s[l], MSL was the first since Viking to use an offset 
center of gravity to fly a lifting trajectory. Like the Viking 
entry vehicles, a reaction control system was used to control 
the vehicle. In addition, MSL was the first to use the RCS 
system to fly a guided entry to land within a small landing 
ellipse on Mars. During the planning for and execution of 
the flight reconstruction, it became obvious that the flight 
mechanics of the MSL vehicle with its offset center of gravity 
are more complex than those of more recent Mars missions 
which had all flown ballistically with no eg offset. The MSL 
data reduction was the inspiration for the analysis in this 
paper and will be used as an example [2], 

The full equations of motion provide the tools to under- 
stand the impact of the capsule flying at an asymmetric 
orientation relative to the approaching wind. The yawing 
moment equation is separated into a linearized system of 
equations for a blunt body flying at angle of attack. The 
re-expression of the moment equation is used to show how 
the yawing moment slope, C np , is properly extracted from 
body rates, r and p, measured by an instrument like an 
onboard inertial measurement unit (IMU). The equations 
show that the natural pitch and roll oscillations are correlated, 
connected by the yaw and roll moment slopes, C np and Cj a 
and the moments of inertia about the roll and pitch axes. A 
flight simulation of the MSL entry is used to validate the 
method of extracting the yaw stability and reconstruction 
results of the MSL flight trajectory are presented as well. 
These cases provide examples of how a simple relation can 
accurately extract the yaw stability for a lifting blunt body 
from measured body rates. 

Like the static moments, the pitch and yaw damping coeffi- 
cients should be modeled asymmetrically for a lifting flight. 
Unfortunately, it is very difficult to measure or compute 
the separate contributions to damping from rotational rates 



d = 4.500 m d = 4.519 m d = 4.520 m 

S = 15.904 m 2 S = 16.039 m 2 S = 16.048 m 2 

Mass: m = 3152.5 kg 

Figure 1 . MSL entry vehicle dimensions 


about the body and wind-relative axes. This will be briefly 
discussed below. For blunt bodies the damping moments 
are small compared to traditional aircraft which have long 
moment arms from the center of gravity and large aerody- 
namic surfaces that resist rotational motion. While small, 
the damping characteristics can have significant impact on 
blunt vehicle dynamics, especially for vehicles without a 
control system. Analysis described briefly below found that 
better measurements or predictions of blunt body damping 
derivatives are required before a correct implementation of 
the damping terms in simulation will have any significant 
influence of capsule damping predictions. 

Mars Science Laboratory Entry Vehicle 

The Mars Science Laboratory entry, descent, and landing 
(EDL) reconstruction was the genesis of the analysis pre- 
sented here. Therefore, the MSL entry vehicle geometry is 
used here as an example of a blunt body of revolution flown 
at angle of attack where some axisymmetric assumptions are 
no longer valid. The MSL entry vehicle consists of a 4.5 
meter, 70° sphere-cone heatshield with a backshell made 
of a series of truncated cones. The vehicle geometry and 
reference parameters are described in Figure 1 and the aero- 
dynamic coefficient and coordinate system definitions are 
shown in Figure 2. The capsule dimensions shown in Figure 
1 are simplified, but sufficient for describing the geometry 
of the flight vehicle. The design and measured diameters 
(and projected areas based on those diameters) show how 
the capsule design and what was fabricated changed slightly 
from the original specifications, set early in the project. The 
mass listed here is for the fully fueled entry vehicle, mea- 
sured before launch. Some representative details of the MSL 
entry are described in later sections as the lateral-directional 
equations of motion are applied to simulation and flight 
data. For details of the vehicle’s aerodynamic characteristics 
and the reaction control system used to control the vehicle 
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Key Geometric Locations 
(in aerodynamic coordinate frame, ac) 


Nominal 

nose — 0.0 
ynose — 0.0 
Z-nose —0.0 


Measured 

x nose —0.0 
ynose —0.0 
Znose — 0.0 


x r Jd = -0.300 
y cg /d = 0.0 
Z cg jd = -0.0215 


x cg /d = -0.295 
y cg /d =0.0 
Z cg /d = -0.0221 


Figure 2. Blunt body coordinate system showing rates about 
body and wind axes. 


and perform guided entry, refer to References [2], [3], [4], 
[5]. For a full report of the trajectory reconstruction and 
comparisons with preflight predictions of the aerodynamics, 
see References [6] and [7]. 

The moments of inertia of the entry capsule are 


Figure 3 shows the angle-of-attack and angle-of-sideslip 
history of a preflight simulation of the MSL entry vehicle. 
Note that the trim angle starts near —16° (Mach > 30) and 
grows in magnitude to near —20°. The vehicle flies in a 
nose-down attitude because the lift of a blunt body comes 
mainly from the axial force component. Pointing nose-down 
directs the axial component of aerodynamic force upwards, 
producing positive lift for guided entry. Near Mach 2.5, 
the vehicle jettisons six balance masses to move the eg to 
the vehicle centerline and null the trim angle in preparation 
for a safe parachute deploy below Mach 2.0. The mass 
jettison event was colloquially referred to as the “straighten 
up and fly right” or SUFR maneuver. The sideslip trim angle 
remained close to zero. 

During entry the guidance system commanded RCS jet fir- 
ings in two primary modes. The first is Range Control 
where the reaction control system performs a series of bank 
reversals to modulate the lift vector to fly to an accurate 
downrange position. Near Mach 6, the vehicle entered the 
Heading Alignment phase where the control logic steered the 
vehicle directly toward the desired landing location, cleaning 
up any crossrange error incurred during the Range Control 
phase. The discrete moment impulses imparted by the 
reaction control jets when called to perform bank reversals 
and damp rates are detected by the IMU which complicates 
the data reduction process. This limited the areas where the 
stability derivatives could be extracted. The details of the 
data reduction process will be described in a later section. 

The aerodynamics of blunt bodies used to fly a lifting trajec- 
tory are somewhere between the non-lifting, ballistic-entry 
variety and a full aircraft or spaceplane like the space shuttle 
Orbiter (as a point of reference, MSL flew at low lift-to-drag 
ratio of 0.24). Because the vehicle geometry is axisymmet- 
ric, the aerodynamic database is still generated as a function 
of total angle of attack (populated with computational fluid 
dynamics codes and ballistic range and wind tunnel data) 
and decomposed into components, but as will be shown 
there is lateral-directional coupling when flying at angle of 
attack much like an aircraft might see. The coupling is 
much more tenuous however, so the equations of motion can 
be simplified to yield practical analysis tools and relations 
between roll and yaw motions which would not be so “clean” 
for an aircraft. While the equations of motion used for the 
analysis done here are straight from classic documents, the 
authors believe this is the first explicit analysis of the static 
stability in this middle ground. 


I xx = 4800 kg — to 2 

I yy = 3800 kg — to 2 (1) 

I zz = 2900 kg — m 2 


1 xyi d-yz ^ 0.01 * I zz 

I xz < 0.04 • I zz 


( 2 ) 


While approximate, these are representative of the values 
used in both the preflight simulation used later for validation 
and the measured values of the flight vehicle. Note that 
all cross products of inertia are small with the largest terms 
being less than four percent of any of the components in the 
trace of the inertia tensor. Most of the terms are less than one 
percent. This permits the longitudinal and lateral-directional 
equations to be decoupled and simplifies the terms of the 
equations developed below. 


3. Equations of Motion 

This section provides a brief overview of how the 6 degree- 
of-freedom (DoF) equations of motion were arranged to 
extract the yaw and roll stability with sideslip angle, /3. 
There are many methods that have been used to identify 
aerodynamic coefficients from flight data. The method for 
extracting static stability described here was selected as 
it uses rate measurements taken directly from an onboard 
inertial measurement unit with no need for further integration 
or manipulation. No explicit attitude information is required, 
yet accurate values of the pitch and yaw stability can be 
measured and compared to preflight predictions to assess 
their accuracy. The rate data are used to extract static stability 
coefficients with very simple relations using only measured 
mass properties and the dimensions of the vehicle with the 
reconstructed dynamic pressure. The method is tolerant 
of many simplifying assumptions and can obtain accurate 
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results very quickly and easily. When a vehicle trims at 
angle of attack, the process is a little more complicated and 
the simplifications introduce small errors in the reconstructed 
yaw stability. However the equations described here are 
useful tools for any blunt body entry data reduction process. 

Assumptions 

The equations of motion used in the derivations shown here 
and below are taken from the work of Duke, Antoniewicz 
and Krambeer [8]. There are a number of simplifying 
assumptions that can be made for blunt bodies that enable 
simple but accurate solutions to the equations of motion for 
use in extracting the stability coefficients. The equations 
shown here are linearized about the trim angle. For non- 
lifting bodies the trim angle is aTrim = P Trim = 0° 
and for a lifting vehicle the equations are linearized about 
a non-zero trim angle, a = aTrim and p Trim = 0°. The 
amplitude of oscillation is assumed to be small. Also, as 
was shown to be valid for MSL, it is also assumed that the 
cross products of inertia are negligible. These assumptions 
allow the longitudinal and lateral-directional equations to be 
decoupled. For a blunt body traveling at hypersonic and 
supersonic speeds, gravity effects on the body rates are small 
over short segments of the trajectory. The change in dynamic 
pressure (or velocity and density) as a vehicle decelerates 
and descends through the atmosphere is assumed small over 
one or two cycles of oscillation. The contributions to angle 
of attack and sideslip due to the transverse accelerations 
from aerodynamic forces, including wind gusts are assumed 
small. The terms which contribute to these accelerations 
will be shown for completeness, before showing solutions 
with those terms omitted. The dynamic stability coefficients 
are assumed to be small, so terms in the moment equations 
proportional to the body rates are neglected. Over most of the 
trajectory the MSL vehicle was damped slightly or saw close 
to neutral dynamic stability and saw very little oscillation 
amplitude growth or decay during entry. Regardless, pitch 
damping does not contribute significantly to the oscillation 
frequency which is of primary interest here. 


Longitudinal and General Non-lifting Capsule Oscillations 

An axisymmetric blunt body with the center of gravity on 
the centerline and forward of the neutral point is statically 
stable and will tend to return to ar = 0° when displaced 
by any total angle-of-attack. As a convenience, an arbitrary 
coordinate system is typically placed on this axisymmetric 
vehicle, defining a pitch and yaw plane with three orthogonal 
forces and three orthogonal moments used to describe the net 
effect of aerodynamic pressures acting on the vehicle. If the 
eg is offset from the centerline, the vehicle will trim at an 
angle away from = 0° and produce lift. Refer back to 
Figure 2 for the coordinate definitions, angle-of-attack and 
angle-of-sideslip definitions, and aerodynamic coefficient 
definitions for the Mars Science Laboratory entry vehicle. 
Note that MSL flew with a radial eg offset. A nonlifting 
vehicle like a Mars Exploration Rover entry capsule has 
identical definitions, but with no radial eg offset. 

For a vehicle with no eg offset, the equations describing the 
pitch and yaw moments are of the same form. For example 
the moment equation about the pitch axis is 


a = q — 


^oc s 

2m 


C La a 


(3) 


Q = 


Qoo Sd 



a 


(4) 


Omitting the heaving term in Equation 3 simplifies this sys- 
tem of equations into a single equation, that of an undamped 
harmonic oscillator. The solution becomes 


a = a 0 cos(u> q t + 5) 


( 5 ) 


The frequency of oscillation is equal to 


4 


(6) 


/9oo Sd 

— \ / t ^rr 


Which can be rearranged to extract the static stability 

_ lyy^q 
l ''m a 0 , 

qocod 

A similar equation can be derived for the yaw stability 


_ I zz u 2 r 

qooSd 



Where ui q and ui r are the frequencies of oscillation of the 
pitch and yaw body rates. For small trim angles and small 
oscillations, these stability relations are of identical form 
and are equivalent. Differences in moments of inertia about 
the two axes would be the only reason the pitch and yaw 
frequencies might differ. For past Mars missions flying non- 
lifting trajectories [9], [10], [11], Equations 7 and 8 provided 
useful tools to assess the static stability and its variation 
along the trajectory using the body rates only. 

For Mars Science Laboratory which flew at a non-zero trim 
angle the solution to the pitch stability equation is nearly 
identical to Equation 7. The difference being the oscillation 
amplitude is centered about the trim angle. The frequency 
of oscillation is still governed by the local pitch stability. 
The yaw equation becomes more complicated and will be 
described in the next subsection. 

Lateral-Directional System of Equations 

At angle of attack, the 3 axis of rotation is no longer aligned 
with the body yaw axis. Figure 4 shows the coordinate 
definition plot again, this time shown at the trim angle 
(3 = 0°) with vectors representing the body and wind- 
relative rotation rates. Note the misalignment between $ and 
the body rates for the capsule at an angle of attack. This 
misalignment introduces a complication for extracting the 
yaw static stability coefficient, C n „. A vehicle like MSL 
typically flies with an onboard IMU and records the body 
rates, p , q , and r. In the case of MSL this data is used for 
real-time navigation during entry and saved for post-flight re- 
construction. Angles of attack and sideslip are reconstructed 
using rate and acceleration data and in the case of MSL, 
surface pressure measurements. Here a system of equations 
describing the yaw and sideslip moments is solved to show 
how the roll, yaw and sideslip rates are all correlated with 
each other. This correlation is used to determine relations for 
identifying the roll and yaw stability coefficients. 

Here a system of equations is defined describing the angular 
accelerations about the roll and yaw axes and the relationship 
between the sideslip and body rates. This matrix was derived 
from the equations described in the NASA document by 
Duke, Antoni ewicz and Krambeer [8]. The equations are lin- 
earized about a trim angle of attack ( otTrim = a 0 , 3Trim = 
0). Note the coefficient subscripts in this matrix equation 
indicate partial derivatives with respect to /3, p, and r. 


Figure 4. Body and wind relative rates 
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The intent here is to develop equations relating the frequency 
of oscillation to the static moments. As mentioned above, 
the damping coefficients are neglected as are the gravity and 
heaving terms. The terms remaining in these equations are 
then 


3p = (sin 3o{C Nf) sin a Q + C A? cos a 0 - C Y ) 

+ cos 3 o(C n sin a 0 + Ca cos a Q + Cy a )) « 0 



(10) 

3 P = sin a 0 

(11) 

3r = — COS a 0 

(12) 

. QooSd 

P0 = r C b 

J-XX 

(13) 

Pp~ 0 

(14) 

Pr ~ 0 

(15) 

. Qoo Sd 

rp = T C nf3 

*zz 

(16) 

fp « 0 

(17) 
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r 


( 18 ) 


Setting the contributions to $ from the aerodynamic forces in 
Equation 10 to zero is an important assumption and required 
to obtain the clean solutions for roll and yaw stability below. 
Linearized about /3 0 = 0°, the first term is clearly zero. How- 
ever, the contributions to /3 from the axial, normal and side 
forces, while small for a blunt body, would provide a more 
accurate expression for sideslip if accounted for. Neglecting 
this term does introduce some error to the expression for 
yaw and roll stability. This will be discussed below when 
comparing data reduction results from a validation case, but 
in short, this simplification does not add a large error. 

These linear, state-space equations for the lateral-directional 
dynamics were transformed to a transfer function represen- 
tation by applying the Laplace transform to Equation 9. The 
transfer function expression was conveniently solved through 
the inverse Laplace transform, yielding a time varying ex- 
pression for sideslip and lateral-directional body rates. The 
resulting solution takes the form of a sinusoid with frequency 
ui, as given by 


Next, the 2nd derivatives of these body rates in Equations 
20 and 21 are equated with the first derivatives of the ex- 
pressions of p and r in Equation 9. Using the expression 
for p in Equation 19 and rearranging yields two equivalent 
equations for the square of the natural frequency of pure 
sideslip oscillation: 


p = -Cip/3 = s f^Ci /3 (p sina 0 -r cosa 0 ) 
= — UJ 2 p + CiU J 1 = W 2 (c 1 — p) 


(25) 


and 


2 qooSd ( psin a 0 — r cos a c 

j I — 

I XX \ ^1 P 


(26) 


r = q f^C n/3 $ = ct f^C np (ps\na 0 - rcosa 0 ) 
= — cu 2 r + c 2 cu 2 = w 2 (c 2 — r) 


(27) 


/ 3(t ) = p o cos cut + 


/ p 0 sin a, 

V w 


r 0 cos a 0 


cu 


sin cut (19) 


ui = 


qooSd 


C, 


np 


p sin a Q — r cos a Q 
c\ — r 


(28) 


P(t) = ^^C uamwt 


p 0 cos Ot- O qoo SdCna r Q cos a 0 qooSdCi a 

Po I J — ^2 


p 0 cos cioQooSdCnp r Q cos OLoqooSdCip 


cos Ult 


(20) 


r{t) = Po i^ d C np sin ut 


Po sinaoqooSdCnp r 0 sin a 0 qooSdCi ^ 


I z Z U 2 


IxxU 2 


cos cut 


Po sinctoqoaSdCnp r 0 sin SdCi^ 


IzzU 2 




( 21 ) 


Where p o , p 0 , and r Q are the initial sideslip, roll and yaw 
rates and lu is the oscillation frequency equal to 


Qoo ^ d 


q<x> Sd 


CJ = \l — C nf) cos a 0 c h sin a Q (22) 

±zz *xx 

The roll and yaw rate equations are of the form: 


Equations 26 and 28 are two expressions for the same 
frequency of oscillation common to the roll and yaw rates. 
Equating the two RHS terms in these equations and rear- 
ranging shows how the static stability coefficients and mass 
properties are correlated with the body rates. 

r-c 1 = Cn IL I (29) 

P - Cl 6 I zz 

The constant terms in Equations 20 and 21 are equal to the 
coefficients c\ and c 2 . Comparing these terms reveals that 
the two constants are proportional to each other. Introducing 
another arbitrary constant, C3 yields 


ci = C3 cos a 0 (30) 

C2 = C3sino! 0 (31) 


Therefore, the constant components of the p and r rates that 
satisfy the system of lateral-directional system of equations 
are correlated. Substituting these rates into the p relation 
within Equation 9 shows that these constant rates are or- 
thogonal to the sideslip rotational vector and represent the 
components of a pure bank (cf>) about the velocity vector. 
Therefore the terms, Ci — p and c 2 — r represent oscillatory 
components of the roll and yaw rates. 


p = a\ cos cut + bi sin cut + Ci (23) 


Pose — P Cl 


(32) 


r = a 2 cos cut + b 2 sin cut + c 2 (24) 


r OS c = r - c 2 


(33) 
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Note that any combination of persistent roll and yaw rates 
that do not meet the orthogonality relationship in Equations 
30 and 31 would result in angle of attack oscillations and 
the longitudinal and lateral-directional equations would no 
longer be decoupled. For a vehicle that is banking, constant 
rates can be subtracted from the measured body rates to 
determine the oscillatory components. Once the constant 
rates have been subtracted, the remaining oscillatory parts of 
the roll and yaw rates are the components of a pure sideslip 
oscillation. Equation 29 becomes 


f osc C'np ^xx 

Pose f '/ Izz 


(34) 


This expression is quite powerful. It says that the roll and 
yaw rates are correlated by the yaw and roll static stability 
coefficients and the moments of inertia about those axes. For 
a sideslip oscillation, the yaw and roll rates are in phase and 
differ by a constant scale factor. The stability coefficients can 
then be extracted from these body rates. For a non-banking 
case. Equation 28 is rearranged into an expression for the 
yaw stability: 


a 


n , 3 


Izz^j T osc 

Qoo Sd 


qooSd (§ 2 ^ sina D - cos a Q ) 

1 osc 

(35) 


and likewise for roll stability 


sideslip rate will alter and complicate this equation, more 
accurately capturing the relationship between the body and 
wind-relative oscillations, but complicating the data reduc- 
tion process. 

Yaw Damping 

In addition to the static moments being complicated by the 
non-zero trim angle of attack, consider how the dynamic sta- 
bility of a blunt body is modeled. The yaw damping moment 
has typically been modeled with the following relation (or 
something similar) for blunt bodies, regardless of trim angle. 


N , 


damping 



c n .) 


(r — /3)d 
4Voo 


(38) 


As it is very difficult to separate the damping due to yaw rate 
from the damping due to sideslip rate, the damping coeffi- 
cients have historically been lumped together and multiplied 
by the average of the yaw and sideslip rates (or one of the two 
selected over the other) to calculate the damping moments. 
As has been shown above, for a vehicle that trims at angle of 
attack, the axes of rotation for these two rates are not parallel, 
so this averaging is an additional approximation. 

This issue has long been a problem for high angle-of-attack 
aircraft dynamic testing and modeling [12], [13]. For exam- 
ple, forced oscillation tests typically impose an oscillating 
yaw rate, r, to a wind tunnel model and extract a component 
of the lumped term 


d xx UJ Pose _ dxx^ 1 

3 q^Sd qooSd (sin a 0 — C os a 0 ) 

(36) 

Analysis of typical blunt capsules shows that the sin a 0 
term is small compared to the cos a 0 term for reasonable 
blunt body trim angles. For a very useful first order approx- 
imation, the yaw stability is scaled by the cosine of the trim 
angle for a given measured frequency. 


C nr =C nr -C n . cos a (39) 


A common practice for aircraft aerodynamic modeling is 
to drop the C nfj term with the assumption that 3 is small 
[12], This is consistent with an aircraft flying a coordinated 
turn where sideslip remains small as a turn is executed. 
Any extreme maneuvers are transient and not oscillatory, so 
even cases that see large sideslip rates do not experience 
oscillations over a long time where damped or undamped 
vehicle characteristics can drive oscillation growth. 


c _ IzzU 2 
3 qooSdcos a 0 


(37) 


This development of equations is shown to provide rigorous 
explanation of how the oscillatory rates are correlated and 
how constant body rates can be present as well. Equation 37 
could also have been determined directly from Equation 22 
by making the assumption that the C) a term was small. Also 
note that Equation 37 reduces to Equation 8 when a 0 = 0°. 

Theoretically, the frequency can be extracted from a roll or 
yaw rate history measured as the frequency of both are the 
same and in phase. In practice, the yaw rate, r , has a greater 
amplitude, thus providing a stronger signal for more accurate 
identification of aerodynamic parameters. Equation 37 is 
a very simple relation between the yaw stability coefficient 
and the oscillation frequency of the body rate data, most 
easily extracted from the yaw rate, r. This equation for 
yaw stability is equal to Equation 8 scaled by the inverse of 
cos a 0 . Accounting for heaving and gravity effects on the 


In contrast, blunt entry vehicles have no large aerodynamic 
surfaces to damp rates. They can see persistent oscillations 
about the pitch and yaw axes that can produce significant 
amplitudes depending on the damping characteristics (often 
dynamically unstable at supersonic speeds). Therefore the 
C' n ,, term can not be neglected especially at low speeds and in 
a windy environment where wind-relative and body rates can 
differ considerably. An extensive study was undertaken to 
improve how the yaw and sideslip damping coefficients were 
modeled for a vehicle like MSL, trimming at angle of attack. 
The analysis was hampered by the long-existing problem of 
having no practical method of separating the damping co- 
efficients. While a mathematical expression was developed 
which more correctly implemented the damping coefficients 
similar to the static relations developed above, in simulation 
it was shown that the new mathematical formulation did not 
affect capsule dynamics appreciably. Until there is a method 
to separate C nr from C n ^ the authors found no compelling 
reason to deviate from the yaw damping implementation 
shown in Equation 38 even when flying at an angle of attack 
as MSL did. 
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4. Data Reduction 

Comparisons to preflight simulation and flight data are de- 
scribed below to show how the modifications to the axisym- 
metric assumptions enable the correct extraction of lateral 
stability of the capsule. Refer to papers by Karlgaard et 
al. and Schoenenberger et al. [6], [7] for a more detailed 
description of the reconstruction and comparisons with the 
preflight aerodynamic predictions. 

Static Stability Identification Process 

For the data reduction process used to compute C ma and 
C np , segments of the vehicle yaw rate data were extracted 
at points along the trajectory where RCS jets were not firing. 
The selection of these quiescent periods ensured that the data 
being used to identify aerodynamic properties of the vehicle 
were taken from regions of the trajectory where the vehicle 
was subject to only aerodynamic forces (and gravity). In 
order to compute C np , a range or window of yaw rate data 
was taken around each of these points and a sinusoid function 
of the following form was fitted to this “truth” signal: 


r fitted = A cos(o ;t + 5) + r c (40) 

In this equation, A is the sinusoid amplitude, u> is the 
frequency of oscillation, <5 is a phase shift and r c is a constant 
offset representing any constant component of the yaw rate. 
A similar equation was fit to the pitch data to compute 
C ma and the angle-of-attack data to obtain the trim angle, 
a 0 , used to compute C np in Equation 37. A simple cost 
function was constructed to take the root sum square of the 
difference between the true yaw rate signal and the fitted 
signal calculated from Equation 40. An optimizer was used 
to find values of A, u>, S and r c that minimized the cost 
function producing the best fit to the true yaw signal. The 
length of the data window over which the signal was fitted 
was chosen to provide one to two oscillations of the sinusoid. 
Through this cost function minimization technique, the fitted 
pitch and yaw frequencies of oscillation were extracted and 
used to compute C' m „ and C n „ . 

Simple Simulation 

Before data from a full entry simulation was used, a simple 
case was run approximating a ballistic range shot. Gravity 
was set to zero and a constant coefficient aerodynamic model 
was used so the pitch and yaw stability would be invariant 
across the simulation. This simulation was intended to 
provide a set of body rates generated from a trajectory in 
agreement with all the simplifying assumptions made in the 
development of the data reduction relations. The extracted 
parameters agreed very well with constant stability coeffi- 
cients in the aerodynamic model. The change in dynamic 
pressure over the length of the sinusoidal wave segment fit 
to the rate data was one notable source of error remaining 
in the data reduction process. However, the error introduced 
by finding the mean frequency over a finite length of time 
was determined to be small. This method of identifying the 
frequency of oscillation from the rate data was selected as a 
reliable and robust technique. There are other methods that 
may produce more accurate results at a discrete point. Such 
improvements are left for later work. 

Validation Simulation 

A representative MSL trajectory was selected to validate 
the relationships between pitch and yaw static stability and 



Figure 5. Mach and dynamic pressure profiles for MSL 12- 
GAL-01 simulation 


the frequencies of oscillation given in Equations 7 (for a 
lifting vehicle) and 37. A preflight simulation provided 
flight dynamics, reaction control inputs and atmosphere and 
gravitational effects that very closely matched the actual 
flight trajectory with the benefit of having perfect knowl- 
edge of how the aerodynamic characteristics were modeled. 
Therefore, the extracted pitch and yaw stability values could 
be compared to “truth.” 

Figure 5 shows the Mach and dynamic pressure profile for 
the preflight simulation 12-GAL -01. The name uses an 
internal MSL project convention and indicates this was the 
first simulation run in 2012 of a landing in Gale crater. Most 
conditions of this simulation are very close to those predicted 
to occur on the day of landing, so analysis of this simulation 
is representative of what was to be expected (and what in fact 
was measured) from flight telemetry. Early in the trajectory, 
Mach number is increasing both because of temperature 
variations in the rarefied segment of the atmosphere and be- 
cause gravity is causing the capsule’s planet-relative velocity 
to increase prior to entering the more dense “continuum” 
portion of the atmosphere. The peak deceleration occurs near 
Mach 16. 

Figure 6 shows the roll, pitch and yaw body rates for the 12- 
GAL-01 simulation. The most obvious features of this plot 
are the bank reversals, indicated by the coordinated spikes 
in the roll and yaw rates. These maneuvers are executed 
during the range control phase of entry to modulate the 
lift vector in order to reach a targeted downrange location 
before entering the heading alignment phase. Starting near 
690 seconds, the capsule enters heading alignment. Note 
the segments of small but constant roll rate, varying in sign 
during this phase. This is indicative of small residual bank 
rates as the capsule drifts between bank angle deadbands. 
RCS corrections fire when the capsule reaches a deadband 
leaving a residual bank rate in the opposite direction which 
is countered when reaching the other bank angle deadband. 
The large bank rate near the end of the entry profile is a 
180° bank reversal performed to place the powered-landing 
descent-stage’s radar system in the proper orientation for 
operation after heatshield separation and later separation of 
the rover/descent-stage from the aeroshell for landing. Also 
note the variation in pitch and yaw rate oscillation frequency 
with time. The variation is driven primarily by the dynamic 
pressure variation. 
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The rates are shown here to show the raw data from which 
the parameters were extracted. As mentioned earlier, RCS 
firings cause large disruptions to the rate signals. The most 
obvious disruptions are the large bank reversals, but any 
small firing to damp rates or correct the bank orientation is 
enough to confuse the data reduction algorithm. Therefore, 
only segments containing a cycle or two free of RCS firings 
were considered for extracting the yaw stability. Nine seg- 
ments of data were identified. 

Figures 7a and 7b show two fits of the sinusoid function 
(Equation 40) to the yaw rates measured during two of 
the nine RCS-free trajectory segments. These fits show 
how closely the model agrees with the simulated rate data. 
Overall fits to each segment were very good and the extracted 
frequencies are very close to the mean frequencies in the 
simulation even when the amplitudes at each peak did not 
match. Based on differences between the body rate and fitted 
curve wavelengths measured at the mean amplitude, errors 
in frequency using this method range from less than 0.5% 
up to near 2% depending on the particular data segment. 
This results in a small error in the extracted parameters. A 
better functional form, perhaps accounting for the change in 
frequency due to deceleration, and/or amplitude change due 
to damping [14], might identify the frequency of oscillation 
more accurately. The constant amplitude sine wave was 
selected for robustness at the expense of the small residual 
fit error. 

Figure 8 shows the extracted pitch and yaw stability coeffi- 
cients identified at the centers of the nine data segments along 
the 12-GAL-01 simulated MSL entry trajectory. The local 
stability coefficients, calculated by finite difference from the 
MSL aerodynamic database, along the entire entry are plot- 
ted and the nine truth points corresponding to the extracted 
points are noted by symbols. Overall the agreement with the 
simulation data is very good. The extracted coefficients are 
within the same error margins for both the pitch and yaw 
stability. Some consecutive yaw stability coefficients ex- 
tracted at points that occurred during the range control phase 
of entry (Mach 14.9, 8.2 and 6.6) were each approximately 
3.5% lower than the database. Subsequent analysis showed 
that omitting the gravity terms in the 0 equation in this region 
resulted in the most noticeable error when compared with the 




Figure 7. Comparison of raw rate data to parameter 
identification model 
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Figure 8. Validation of static stability method using MSL 
12-GAL -01 simulation 
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simulation (3 history. The heaving terms did not appear to 
appreciably alter 0 for the MSL vehicle. Keeping the gravity 
terms would produce a more accurate solution to the system 
of equations, but the extraction of the yaw stability would 
no longer be dependent on the oscillation frequency only. 
A more complex analysis would be required to extract the 
yaw stability. The errors in this region provide an estimate of 
the accuracy limits when using the simplified relation. The 
correction to Equation 8 reflected in Equation 37 is the first- 
order effect of an axisymmetric capsule flying at a non-zero 
trim angle. 

Mars Science Laboratory Flight Data 

In the same manner as the validation case, the static stability 
coefficients were extracted from MSL flight data using Equa- 
tions 7 and 37. Figure 9 shows the results as published in a 
paper by Schoenenberger et al. [7]. The RCS firing history 
during MSL’s entry was less active than preflight predictions 
over much of the trajectory. However, there were small but 
frequent firing events early in the trajectory which meant 
fewer quiescent data segments in the hypersonic regime. 
Figure 10 shows representative examples of the rate data 
measured by the Descent Inertial Measurement Unit (DIMU) 
during entry compared to the data reduction functions to 
which the data were fit. The rate data were very clean and 
fitting the sinusoid was much like doing so with simulation 
data. 

The data points in Figure 9 are plotted against stability 
coefficients predicted by the MSL aerodynamic database, 
queried along the reconstructed trajectory determined by 
Karlgaard et al. [6]. Agreement with the extracted data and 
the nominal database is excellent with errors very similar to 
the validation case. This indicates that flight data quality was 
very good and the preflight database closely predicted the 
stability characteristics of the capsule during all of entry. 

As MSL trimmed at an angle of between aTrim = —16° 
and —20°, the correction to the axisymmetric equation for 
C nj3 removed an error of between 4 and 7% across the entry. 
This is important when evaluating the reconstructed data 
and comparing to preflight predictions and the uncertainty 
margins used in preflight trajectory and controller design. 
This improvement removes an error that would have been 
a significant fraction the of the ±20% static stability uncer- 
tainty carried in the MSL aerodynamic database [4]. Without 
the understanding provided by the formulations shown here, 
the MSL project may have recommended expanding the yaw 
stability uncertainties. 


5. Conclusions 

The Mars Science Laboratory entry vehicle was the first 
planetary spacecraft to fly a lifting trajectory since Viking 
and the first to fly a guided entry at another planet. Prior 
to MSL a series of axisymmetric vehicles flying ballistic 
trajectories were flown to Mars. Assessment of the flight 
performance of these vehicles was simplified because of the 
axisymmetric shape of the vehicle and the trim angle being 
close to a total angle of attack of 0°. The MSL capsule 
trimmed between —16° and —20° during its entry through 
the Mars atmosphere. Flying at this angle effectively changes 
the aerodynamic characteristics into those of a vehicle some- 
where between an axisymmetric capsule and a traditional 
aircraft like the space shuttle Orbiter. One key aspect that 
changed was the static stability. Instead of describing the 


Mach (approximate) 



Velocity (m/s) 


Figure 9. Comparison of flight stability coefficients to 
aerodatabase [7] 



Figure 10. Comparison of parameter ID fit to measured 
MSL body rates at Mach 3.6 [7] 


blunt body aerodynamics with two separate and orthogonal 
stability equations of the same form, MSL required longitu- 
dinal equations to describe the pitching motion and lateral- 
directional equations where roll and yaw motions are cou- 
pled. These equations were solved and several simplifying 
assumptions were made to yield simple expressions for the 
roll and yaw stability in terms of the frequency of oscillation 
common to the roll, yaw and sideslip rates. 

The yaw stability equation was used to extract coefficients 
from a simulated MSL trajectory. Agreement between the 
extracted parameters and the aerodynamic database was 
good, with the omission of gravity terms being identified as 
the most significant source of error in this application. A 
solution of the full set of equations could be done to account 
for the dropped terms at the expense of the simple relation 
between static stability and frequency. Results for the flight 
data reconstruction, first reported in another paper [7], show 
extracted pitch and yaw stability coefficients in agreement 
with preflight predictions to within an accuracy very similar 
to the validation case run for this paper. The aerodatabase 
predicted the MSL static stability very well. 

A brief discussion was included regarding how damping 
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coefficients might be more accurately implemented in sim- 
ulations. As the symmetric static stability of blunt bodies of 
revolution is broken when flying at angles of attack, so too 
are the dynamic damping characteristics. This has long been 
an area of research for high angle-of-attack aircraft flight 
and is often limited by the lack of experimental capabilities 
for separating damping derivatives due to body and wind 
relative attitude rates. Analysis done in support of this paper 
indicated that without much better information about the 
differences between these damping derivatives, there is no 
significant benefit to altering the combined coefficients as is 
done with axisymmetric vehicles. 
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